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Abstract 

Commonly studied cellular automata are memoryless and have fixed topology of connections be- 
tween cells. However by allowing updates of links and short-term memory in cells we may potentially 
discover novel complex regimes of spatio-temporal dynamics. Moreover by adding memory and dy- 
namical topology to state update rules we somehow forge elementary but non-traditional models of 
neurons networks (aka neuron layers in frontal parts). In present paper we demonstrate how this 
can be done on a self-inhibitory excitable cellular automata. These automata imitate a phenomenon 
of inhibition caused by hight-strength stimulus: a resting cell excites if there are one or two excited 
neighbors, the cell remains resting otherwise. We modify the automaton by allowing cells to have 
few-steps memories, and make links between neighboring cells removed or generated depending on 
states of the cells. 

Keywords: cellular automata, dynamical topology, memory, excitation, neural networks 

1 Cellular Automata and their excitable species 

Cellular Automata (CA) are discrete, spatially explicit extended dynamic systems. A CA system is 
composed of adjacent cells or sites arranged as a regular lattice, which evolves in discrete time steps. 
Each cell is characterized by an internal state whose value belongs to a finite set. The updating 
of these states is made simultaneously according to a common local transition rule involving only a 
neighborhood of each cell |14|.|31]. Thus, if u| is taken to denote the value of cell i at time step T, 

the site values evolve by iteration of the mapping: er£ T+1 ' = 4>{a^p £ Mi \ , where cj> is an arbitrary 
function which specifies the cellular automaton rule operating on the neighborhood M of the cell i. 

In totalistic rules the value of a site depends only on the sum of the values of its neighbors, and 
not on their individual values: cr' T+1 ^ = (f>(^ "j^') • 

This paper deals with two-dimensional semi-totalistic excitable CA, where every cells take three 
states: resting 0, excited 1 and refractory 2. State transitions from excited to refractory and from 
refractory to resting are unconditional, they take place independently on cell's neighborhood state: 
o-P = 1 -> a? +1) = 2, a f> = 2 -+ = 0. 

In representing excitation rule we adopted Pavlovian phenomenon of defensive inhibition: when 
strength of stimulus applied to some parts of nervous system exceeds certain limit the system 'shuts 
down', this can be naively interpreted as an inbuilt protection of energy loss and exhaustion. To 
simulate the phenomenon of defensive inhibition we adopt interval excitation rules, developed in [3], 
and put as a resting cell becomes excited only if one or two of its neighbors are excited. If more then 
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Figure 1: Configurations of defensive inhibition rule CA developed from an excited singleton. Resting 
cells are white, excited cells are red and refractory cells are blue colored. 

two neighbors excited the defensive inhibition comes into action and prevents the cell from excitation: 
af>=0^f> = lif £>f> = 1)6 {1,2}. 

An example of simple development is shown in Fig. Q] where configurations of defensive-inhibition 
CA starting from an excited singleton, and evolving up to T — 15; there, as well as Sect. 2 and 
initially in the remaining, the Moore neighborhood (eight neighbors plus central cell) is adopted. 

The defensive inhibition rule was studied in [I] amongst other interval excitation rules in a con- 
text of morphological and dynamical complexities. We demonstrate there that the rule exhibits quite 
a low morphological complexity — it is amongst rules with lowest number of different cluster (i.e. 
clusters of excited and refractory states) sizes. However CA governed by the rule exhibits longest 
transient periods. Calculating Langton's A parameter we also shown that the rule stays toward the 
middle of the rule phases space, and occupies general position of rules with complex behavior. A 
contradiction between morphological, i.e. a posteriori, and function-based A, i.e. a priori, measure- 
ments of complexity sounds alarming. Moreover, the closest to the defensive inhibition rule - in the 
interval excitation universe - is the interval [2,2] rule is a 'kingdom of complexity and universality', 
as well demonstrated in |3]. Mobile localizations, in other terminologies called wave-fragments and 
gliders, are essential attributed of complexity. As you can see on Fig. [1] the rule is a distant anolog 
of replicator rules. 
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Figure 2: Snapshots of 300 x 300-ccll automaton governed by defensive inhibition rule, (a) development 
from excited singleton, configuration at step T = 119, (b) development from random configuration of 
uniformly distributed states, sites in the centered sub-lattice of 200 x 200 cells were assigned values at 
random, configurastion recorded at step T = 72. 

A singular initial excitation leads for formation of four mobile localizations (wave-fragments) trav- 
eling South- West, South- West, North- West and North-East (Fig. [Tjle). These localizations however 
split onto two localizations each (Fig. [Tf) merge together (Fig. [Tg) but then annihilate almost ev- 
erywhere but in concavities of the growing pattern (Fig. [TJi) . The remained domains of excitation 
produce new mobile localization and the process of replication and annihilation continues. 

Eventually, spatio-temporal dynamics becomes so complicated that large-scale patterns emerged 
from excited singleton look - by naked eye - as disordered as configuration developed from initial 
random configuration (Fig. [2}. These patters are comprised of frequently born mobile localizations, 
gliders, which collide with each other and produce new localization. The frequency of their reproduc- 
tion is however so high and the localization's life-time is so short that it is impossible to trace them 
visually. 

2 Cellular Automata with Memory 

Standard CA are ahistoric (memoryless) : the transition function depends on the cells neighborhood 
configuration only at the preceding time step. Historic memory can be embedded in CA by presenting 
every cell by a mapping of its states in the previous time steps. Namely, we propose to maintain the 
transition rules (<f>) unaltered, but to make them acting on the cells featured by a function of their 

previous states: f | T+1 ^ = ^( s j T ' NjJ, where s[ T ' is a state function of the series of states of the 
cell i up to time-step T. Cellular automata implementing memory in cells will be termed historic, 
and the standard ones ahistoric. 

Thus, cells can be featured by the most frequent of their three last states |10j : 

sT^modeiar^.of^^}. 
In case of a tie, i.e. the three states present once, the last state will feature the cell. The historic and 
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Figure 3: Effect of mode memory starting from a simple configuration. The last series of evolving patterns 
shows the underlying patterns with cells featured by their most frequent state (s) along the last three 
time steps; these patterns generate the actual patterns (a) with mode memory above-right them. 
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ahistoric patterns coincide up to T 
becomes operative after T = 3 . 

Figure [3] shows the effect of mode memory starting from a simple configuration. As a general 
rule, memory tends to restrain the evolution as shown in the case of mode memory in Fig[3] It is 
generally so from the beginning of the effective memory action, so at T = 3 the outer excited cells 
in the actual pattern evolution are not featured as excited but as resting cells, as this is the their 
most frequent state up to this time step (twice resting versus one excited). Typically, the series of 
evolving patterns with memory diverges from the ahistoric evolution already at T = 4. From this 
early time-step, the patterns with memory turn out less expanded, as shown in FigfjJ] The restraing 
effect of memory may result fatal starting with simple configurations such as the singleton of FigFJ] 
which extinguishes at T = 4 with mode memory , a fact that led us to change the simple initital 
configuration to demostration purposes to that of Fig. [3] 

Figure [4] shows the effect of memory when initially every cell of the 100 x 100 cell sub-lattice of 
300 x 300 cell lattice receive resting state with probability 0.8 and excited and refractory states with 
probability 0.1 each. The growing excitation pattern is 'guided' by excited wave-fragments. There 
two immediately visible phenomenological differences between spatio-temporal dynamics of defensive 
inhibition CA with memory and without memory. First, excitation-inhibition pattern in CA with 
memory spread two times slower compared to that of CA with memory. So, memory 'delays' space- 
time developments. Second, configurations of CA without memory feature singletons of excited and 
refractory states, while in the configurations of CA with memory, more accentuated fiber-like clusters 
of excited and refractory states. There are also hints that small domains of labyrinthin-like patters 
of refractory states are surrounded by broken wave fronts of excitation (Fig. [4j ■ 

It should be emphasized that the memory mechanism considered here is different from that of 
other CA with memory reported previously (see[TS],p.ll8; [U,p.43; or class CAM [2], p. 7). Typically, 
higher-order-in-time rules incorporate memory into the transition rule, determining the configuration 
at time T + 1 in terms of the configurations at previous time-steps. Thus, in second order in time 
(memory of capacity two) rules, the transition rule operates as : a 



(T+l) 
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1 The only cells not featured as resting after T = 3 are those refractory at T 
extinction. 



which implyes immediate 
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(d) T = 120 



(e) T = 150 
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Figure 4: Development of CA with memory from random small configuration. 
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. Double memory (in transition rule and in cells) can be implemented as : a\ = &[Sj G 
A/i,s^ T 1 ' G A/ij. Particularly interesting is the reversible formulation: <r| T+1 ' = <^o"j T ' £ A/iJ 
ff^ ; , reversed as of ; = Ma) ' G M) G of • To 

preserve reversibility, the reversible 

formulation with memory must be : <r' T+1 ' = <f>(s^ G A/iJ cr^ T ' — 1 ^ [11] . 

Some authors, for example Wolf-Gladrow [17], define rules with memory as those with dependence 
in <j> on the state of the cell to be updated. So elementary rules with no memory, such as rule 90, 
take the form : a\ ' = 4>[ of a\ +1 J . Our use of the term memory is not this. 

CA with memory in cells are cited by Wuensche and Lesser, [18],P-15, who do not enter into its 
study and state that "CA with memory in cells would result in a qualitatively different behavior " . 

Let us point here that the mechanism of implementation of memory adopted in this work, keeping 
the transition rule unaltered but applying it to a function of previous states, can be adopted in any 
dynamical system. In an earlier work f[9].[8].[T0]) we explored the effect of embedding this kind of 
memory into discrete dynamical systems : xt+i = f(xr) by means of xt+i = /(rar) with mj being 
a mean value of past states. We have studied this approach in, perhaps, the canonical example: the 
logistic map, which becomes with memory xt+i = mr + Xmril — bit). In [10], we studied the 
effect of memory in a particular Markovian stochastic processes (the random walk), p' T+1 = PyM 
by means of p^ +1 = 7TyM with ttt being a weighted mean of the probability distributions up to T. 



3 Structurally Dynamic Cellular Automata 

Structurally dynamic cellular automata (SDCA) were invented by Ilachinski and Halpern [19]. They 
suggested the connections between the cells, are allowed to change according to rules similar in nature 
to the state transition rules associated with the conventional CA. This mean that given certain 
conditions, specified by the link transition rules, links between rules may be created and destroyed. 
The neighborhood of each cell is now becoming dynamic rather than fixed throughout the automaton. 
Thus, state and link configurations of an SDCA are both dynamic and are continually altering each 
other. 

In the Ilachinski and Halpern model, a SDCA consists of a finite set of binary-valued cells 
numbered 1 to N whose connectivity is specified by a NxN connectivity matrix in which Xij 
equals 1 if cells i and j are connected; otherwise. Thus, we have A/j = {j/X^J = 1} and 
^(t+i) _ ^0.( T ) g A/f"'Y The geodesic distance between two cells i and j, Sij, is defined as the 
number of links in the shortest path between i and j. We say that i and j are direct neighbors if 
Sij = 1, and that i and j are next-nearest neighbors if Sij = 2, so A/"A/"' T ' = {j j S^p = 2}. There 
are two types of link transition functions in SDCA: couplers and decouplers, the former add new 
links, the later remove links. The set of coupler and decoupler determines the link transition rule: 

Instead of introducing the formalism of the SDCA, we deal here with just an example, in which 
the decoupler rule removes all links connected to cells in which both values are at refractory state 
(Ay = 1 —¥ A'J +1 ' = iff er ] - T ' = 2 and and the coupler rule adds links between all 

next-nearest neighbors sites in which both values are excited (Ay ' = ->■ X)] = 1 iff of ; = 1 

and <rj T ' =1 and j G A/"A/"' T '). In the SDCA treated here, the transition rule for cell states is that 
of the generalized defensive inhibition rule, excitation interval [1,2]: resting cell is excited if a ratio 
of excited and connected to the cell neighbors to total number of connected neighbors lies in the the 
interval [1/8,2/8]. 

The initial scenario of Fig[5]is that of Fig[3]with the wiring network revealed, that of an Euclidean 
lattice with eight neighbor^. No decoupling is verified at the first iteration in Fig[S] but the excited 



2 Moore neighborhood in CA terminology, in which, the generic cell ■ has sixteen next-nearest neighbors : I ■ 
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cells generate new connections, most of them lost, together with some of the initial ones, at T = 3. 
The excited cells at T = 3 generate a crown of new connections at T = 4. 




Figure 5: The structurally dynamic cellular automaton described in Section 3, up to T = 4. 



4 A Structurally Dynamic Cellular Automaton with Mem- 
ory 

Memory can be embedded in links in a similar manner as in state values, so the link between any 
two cells is featured by a mapping of its previous link values: 1^ = Z(AJ - , . . . , Ay ). The distance 
between two cells in a historic model (dy), is defined in terms of the I instead of the A values, so that 

i and j are direct neighbors if dy = 1, and are next-nearest neighbors if dij = 2; N± = {j/d\^ = 1}, 
(t} i m 

and NN\ = {j j d\- = 2} . In our approach here, the memory rule for links (I) is the same 
that of state values (s). Generalizing the approach to embedded memory introduced in Sect. 2, 
the unchanged transition rules (<f) and 4>) operate now on the featured link and cell state values: 

*r- 1} = *(<r ej*) ,^ = tp(i%\ s ?\s<p) mum- 

Figure [6] shows the effect of mode memory in the scenario of Fig[5] Again, memory has an inertial 
effect. Also in what concerns to the link dynamics: the initial wiring network tends to be restored as 
made apparent in the underlying patterns of Fig(5] But inertial effect does not mean full restrain, as 
FigEl shows. This figure, the ahistoric and mode memory patterns at T — 20, makes again apparent 
the preserving effect of memory . 

Figure [S] shows the effect of mode memory at T = 20 starting at random in what concerns to cell 
state levels in an initially eight-neighbours lattice of size 45 x 45 with periodic boundary conditions. 
The wiring of the border cells are not shown in Fig[5]to avoid that the connections of the border cells 
(linked to their opposite ones in the lattice) do mask the whole pattern. The ahistoric pattern shows 
a low density of non resting cells, together with a irregular tessellation which clearly differs from the 
initial wiring network, but with a notable absence of links traversing the pattern. Opposite, the 
historic pattern exhibits a higher presence of excited and refractory cells as well as higher density of 
links, somehow more reminiscent of the initial structure but much traversed by connections creeated 
during the evolution with memory. 

Figure shows the evolution up to T = 50 of the excited cell density, the average number of 
nearest neighbors and next-nearest neighbors per site and the effective dimension: the average ratio 
of the number of next-nearest to nearest neighbors per site (a discrete analogue to the continuous 
Hausdorff dimension), in the simulations of Fig. [5] The evolution in the mode memory model shown 
in Fig.[9]is that of the featured patterns in all the parameters, and also that of the excited cell density 
and nearest neighbors of the actual patterns, corresponding to the no dot marked curves. Again the 

3 Starting from an Euclidean lattice with eight neighbours, the effective dimension is 16/8=2. 
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(7, A 



Figure 6: The structurally dynamic cellular automaton with mode memory. Evolution from T = 3 up to 
T = 6 starting as in Fig[5] The first row of evolving patterns applies to the actual patterns, the second 
row shows the evolving patterns of the featuring states of cell and links. 

inertial effect of memory is shown in the evolution of excited cell density and neighbors averages: the 
ahistoric curves (in red) stabilize in levels more distant (lower) to the initial ones than those of the 
ahistoric. The excited cell density evolution shows a fairly consistent parallelism between the actual 
(dotted) and featured patterns, whereas the actual and featured average number of nearest neighbors 
tend to approach, after an initial perturbation of the former. The general validity of the form of the 
curves shown in Fig. [9] as well as the aspect of patterns at T = 20 has been assessed by running ten 
more simulations with different random initial configurations. 



5 Other memories 



Memory can be implemented only in the cell state dynamics but not in the link dynamics, or viceversa. 
Figure \W\ shows the effect of such partial memory implementations starting as in Fig[5] Memory in 
cells appears much more determinant in preserving the initial features than memory in links. 

An intermediate way of implementing updating will be that of adopting different time scales: cell 
states are updated every time step, while cell links are updated every 10th or even every 100th time 
steps. This will be an analogue of fast and slow variables in non-linear equations. The idea may be 
translated to the way in which memory is implemented, which is planned for a future research work. 

Memory of previous configurations can be embedded in the CA dynamics by means of a weighted 
average of previous states. Thus, a weighted average memory of all previous states can be imple- 
mented as: 



= < 
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Figure 7: The cellular automaton starting as in FigJS]at T=20, with no memory (left) and mode memory 
in both cell states and links. 



and 4 T) = < 



*iV */ 4 T) = °- 5 , after - (T) 



•ij -j -i — - , uiu^i m^ 







if a\j> <0.5 ! + 



The choice of the memory factor a simulates the long-term or remnant memory effect: the limit 

case a = 1 corresponds memory with equally weighted records (full memory model, mode if k = 2), 

whereas a << 1 intensifies the contribution of the most recent states and diminishes the contribution 

of the past ones (short type memory). The choice a — leads to the ahistoric model. Coding states 

1 

as 0, 1, . . . , k — 1, a-memory is effective if a > 



2(k - 1) 



T 



As an example, Fig llll shows the effect of full memory m' T ' = — cr^f 1 in the initial scenario 

t=i 

of FigO Again, memory restrains the evolution of the defensive inhibition rule, but in a rather 
different way compared to the mode memory in Fig[3] Embedding weighted memory typically has 
effect already at T = 2, as in Fig llll bv featuring the initially refractory cells (code 2) as excited (code 

2 + 

1) as being resting (code 0) at T = 2 : 1 = — - — . This kind of, let us say, ad hoc featuring when 
working with more than three states is an inherent consequence of averaging, that tends to bias the 
result to the mean value. This led us to focus on a much more fair memory mechanism: the mode. 

In Fig |12l memory is implemented in an alternative way: cells are exciting if were excited at T 
and/or at T — 1. This implies a divergence in the evolving patterns from T = 3 (as shown comparing 
Fig |12l to Fig. 0, with a higher chance of resting cells to be in contact with excited cells. 



6 Discussion 

In present paper we undertook scoping studies of the effect of memory embedded in cells and links of 
specific structurally dynamic cellular automaton. We shown, that as a rule, memory notably alters 
the ahistoric dynamics. A complete analysis of the effect of memory on structurally dynamic cellular 
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Figure 8: The structurally dynamic cellular automaton starting at random in a 45 x 45 lattice. Patterns 
at T=50 with no memory (left), and mode memory. 

automata (SDCA) is left for future work which will develop a phenomenology of SDCA with memory, 
i.e. the full analysis of the rule space based on the morphological classification of patterns formed, 
the intrinsic parameters (e.g. Langton's lambda, Wuensche's Z parameter), the structure of global 
transition graphs (this would be feasible only in the ID case), the entropy and other dynamics-related 
issues. Potential fractal features are also to come under scrutiny [24] , 

Some critics may argue that memory is not in the realm of CA (or even of dynamic systems), 
but we believe that the subject is worth studying. At least CA with memory can be considered 
as a promising extension of the basic paradigm. A major impediment to modeling with CA stems 
from the difficulty of utilizing the CA complex behavior to exhibit a particular behavior or perform a 
particular function: embedding memory in states and links broadens the spectrum of CA as a tool for 
modeling. It is likely that in some contexts, a transition rule with memory could match the "correct" 
behavior of the CA system of a given complex system (physical, biological, social and so on). 

To substantiate this our initiative we tested how incorporating of the memory will influence 
dynamics of a system, more close to real life. We have to model Belousov-Zhabotinsky (BZ) reaction, 
a classical example of excitable chemical system. Namely, we have chose to simulate BZ reaction in 
sub-excitable mode: where local perturbations do not lead in general to target or spiral waves but to 
compact localized wave-fragments, or mobile localizations, propagating in the reaction space |25] . 

We simulate the excitable chemical system using the two- variable Oregonator model |15].|29| 
modified to describe the light-sensitive analogue of the BZ reaction under applied illumination [121 1 13J : 

— = -{u-u -(fv + 4>) — ■ — ) + D u X7u 
at e u + q 

dv 

where variables u and v represent local concentrations of bromous acid HB1O2 and the oxidized form 
of the catalyst ruthenium Ru(III) respectively, e sets up a ratio of time scale of variables u and v, 
q is a scaling parameter depending on reaction rates, / is a stoichiometric coefficient. Variable (f> is 
a light-induced bromide production rate proportional to the intensity of illumination and controls 
the excitability of the medium. The excitability is modified by light because the excited ruthenium 
catalyst reacts with bromo-malonic acid to produce bromide, which is an inhibitor of autocatalysis. 
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10 20 30 40 50 10 20 30 40 50 

Time Time 

Figure 9: Evolution ol excited cell density, average number of nearest neighbours and next-nearest neigh- 
bours per site, and effective dimension in the simulation of FigJ5] Ahistoric simulation in red, mode 
memory one in blue. 

The activator diffusion coefficient is D u , whilst the diffusion term for v is absent as it is assumed that 
the catalyst is immobilized in a thin-layer of gel. To integrate the system we used a Euler method 
with a five-node Laplacian operator, time step At— 1CP' ! and grid point spacing A x = 0.25, with 
the following parameters: <j>o = 0.076044, e = 0.02272, / = 1.4, q = 0.002. These parameters roughly 
correspond to the region of higher excitability of the sub-excitable regime [25] . 

The waves were initiated by locally disturbing the initial concentration of species u, e.g. a few grid 
sites in a chain are each given a value u — 1.0; see configuration of initial perturbation in Fig. 113b .. 
In memoryless model such a segment-perturbation generates two localized wave fragments each, the 
wave-fragments travel a finite distance preserving their shape like quasi-particles or dissipative solitons 

(Fig.rnb). 

Analyzing various types of memory in the model we found that not only the memory function can 
boost or inhibit spreading wave-patterns but also distort topology of wave- activity. The following 
modes were considered: 

• model Mi: u l = mai(u l , u t_1 , u*~ 2 ) and v* = max{v l , u t_1 , i/~ 2 ), 

• model M%: variable v is not affected but u* = maxiu 1 , u* , u' -2 ), 

• model M3: variable u is not affected but u* = max{v l , 

• model M 4 : u* = (u* + u^ 1 + u'" 2 )/3. and v* = (v* + v^ 1 + d'~ 2 )/3., 

• model M5: variable u is not affected but v t = (i/ + v t ~ 1 + v t ~ 2 )/S.. 

In the model Mi wave-fragments, generated by initial perturbations start to expand. The same 
stimulus in the model M2 lead to different outcomes, one local perturbation leads to formation of 
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Figure 10: The cellular automaton starting with memory only in links (left), and only in cell states. 
Patterns at T = 20 starting as in Fig|5] 
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Figure 11: The defensive inhibition rules with arithmetic mean memory of all previous states. Initial 
scenario of FigEJ The first row of evolving patterns applies to the actual patterns, the second row shows 
the underlying patterns with cells featured by the rounded average state value of all its previous states. 




■ * o 8 [■] 





Figure 12: The defensive inhibition rule featuring cells as exciting if excited at T and/or at T — 1. Initial 
scenario of Fig. 03 The first row of evolving patterns applies to the actual patterns, the second row shows 
the cells as exciting if were excited at any of the two last time steps. 
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Figure 13: (a) configuration of initial excitation-stimulation, (b) localized wave-fragments in BZ model 
without memory, (c) model Mi, (d) model M 2 , (e) model M 3 , (f) model M 4 , (g) model M 5 . 

target wave, while another perturbation is transformed into two expanding but still localized wave- 
fragments. All activity diminishes in model M3. Both perturbations lead to target waves in models 
M 4 and M 5 . 

The SDCA seems to be particularly appropriate for modeling the human brain function, -updating 
links between cells imitates variation of synaptic connections between neurons represented by the cells 
-, in which the relevant role of memory is apparent. Models similar to SDCA have been adopted to 
build a dynamical network approach to quantum space-time physics [2H] . Reversibility is an important 
issue at such a fundamental physics level; a generalization of the Fredkin's reversible construction 
is feasible in the SDCA scenario, which can be endowed with memory as: f' T+1 ^ = 4>(^s^ 6 

A^ (T> ) ea\ T ~ 1) , A<J +1) = ^(^J ) ,sf ) ,4 T) ) Q \-J~ X) G3- Technical applications of SDCA (and 
graph theory) may also be traced, so [27] and [28] . 

Anyway, besides their potential applications, SDCA with memory have an aesthetic and mathe- 
matical interest on their own. The study of the effect of memory on CA has been rather neglected 
and there have been only limited investigations of SDCA since its introduction in the late eighties 
0. Nevertheless, it seems plausible that further study on SDCA (and Lattice Gas Automata with 
dynamical geometry [30]) with memory Q should turn out to be profitable. 

4 To the best of our knowledge, the relevant references on SDCA are |20|-|23|. together with a review chapter in |14| and 
a section in [2] . 

5 Not only in the basic paradigm scenario, but also in SDCA with random but value dependent rule transitions (which 
relates SDCA to Kauffman networks 20 ), and/or in SDCA with the extensions considered in 1231 . such as unidirectional 
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